Gate defined wires in HgTe quantum wells: from Majorana fermions to spintronics 
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We introduce a promising new platform for Majorana zero-modes and various spintronics appli- 
cations based on gate-defined wires in HgTe quantum wells. Due to the Dirac-like band structure 
for HgTe the physics of such systems differs markedly from that of conventional quantum wires. 
Most strikingly, we show that the subband parameters for gate-defined HgTe wires exhibit exquisite 
tunability: modest gate voltage variation allows one to modulate the Rashba spin-orbit energies 
from zero up to ~ 30K, and the effective g-factors from zero up to giant values exceeding 600. The 
large achievable spin-orbit coupling and g-factors together allow one to access Majorana modes in 
this setting at exceptionally low magnetic fields while maintaining robustness against disorder. As 
an additional benefit, gate-defined wires (in HgTe or other settings) should greatly facilitate the 
fabrication of networks for refined transport experiments used to detect Majoranas, as well as the 
realization of non-Abelian statistics and quantum information devices. 

PACS numbers: 



I. INTRODUCTION 

The ability to efficiently manipulate electron spins 
with electric and magnetic fields underlies a wide vari- 
ety of solid-state applications P Prominent classic exam- 
ples include giant magnetoresistancep^ spin qubitsP® 
and spin transistors! 1 * 7 ! 8 * Recent proposals for stabiliz- 
ing Majorana zero-modes in topological insulatoJ^^ 
and semiconductoi'^K^ architectures, while not usually 
viewed from a spintronics lens, similarly rely crucially on 
spin manipulation. In essence these approaches utilize 
spin-orbit coupling and Zeeman fields to effectively con- 
vert an ordinary s- wave superconductor into a 'spinless' 
topological superconductor supporting Majorana zero- 
modes (for recent reviews, see Refs. 16-19). Intense ex- 
perimental efforts, dri ven p artly by potential quantum 
computing applications^^, have already delivered pos- 
sible Majorana signatures 

For many such spin-based applications, materials ex- 
hibiting easily tunable spin-orbit coupling and g-factors 
are highly desirable. In this paper we employ comple- 
mentary analytical and numerical methods to demon- 
strate that gate-defined wires in HgTe quantum wells 
(see the geometries in Fig. [T]) satisfy both criteria. By 
itself this observation is unremarkable; for instance, 
Rashba coupling in semiconductors is well-known to be 
gate-tunable ™ while g-factors can be modified through 
various means including electric fields and strainP^HsSl 
Rather, the special feature of the HgTe wires we study — 
which stems largely from the unusual Dirac-like band 
structure exhibited by the quantum well — lies in the ex- 
traordinary degree to which these parameters can be con- 
trollably varied under realistic conditions. 

As in any semiconductor, gate voltages can induce 
moderate changes in Rashba coupling for the two- 
dimensional HgTe quantum well hosting the wire. We 
show, however, that the effective Rashba parameters 



for quasi-one-dimensional confined subbands vary much 
more dramatically and in an oscillatory fashion, similar 
to Refs. 31-34, Relatively modest gate voltages can con- 
sequently alter the characteristic spin-orbit energies for 
the wire from zero to appreciable values of ~ 30K (for 
comparison typical spin-orbit energies for electron-doped 
wires such as InAs or InSb are ~ IK 35 ). More surprising 
is the behavior of the effective g-factors for confined sub- 
bands, which in contrast to typical wires are by far domi- 
nated by orbital contributions from the magnetic field (at 
least when directed normal to the well). These g-factors 
similarly undergo gate-induced oscillations and can be 
driven from zero to enormous values exceeding 600 due 
to orbital enhancement. In both cases the remarkable 
oscillatory dependence originates from non-perturbative 
modifications of confined wavefunctions in response to 
gating. 

Because of this exquisite tunability, gate-defined HgTe 
wires are prime candidates for spintronics and related 
applications. Here we focus on one particularly enticing 
example — the pursuit of Majorana modes for topologi- 
cal quantum information processing. (Note that edge 
states of HgTe in the two-dimensional topological insula- 
tor phase can also host Majoranas! 10 * 36 " ^ The physics we 
discuss here is unrelated to these edge states, but is in- 
stead close in spirit to the semiconductor wire proposals 
from Refs. I12|131 ) We show that when a good proximity 
effect with an s-wave superconductor is generated, the 
giant g-factors allow for exceptionally weak fields — a few 
tens of mT — to drive the wire into a topological super- 
conductor with Majorana zero-modes. The strong spin- 
orbit coupling for the HgTe wire (compared to typical 
electron-doped wires) further allows this topological state 
to possess a relatively large gap that exhibits enhanced 
immunity against disorder. 41 Apart from these virtues 
we expect that gate-defined wires offer another impor- 
tant longer-term advantage as well. Namely, synthesizing 



arbitrary wire networks merely requires patterning of ad- 
ditional gates on the quantum well. These can serve the 
dual purpose of enabling refined multi-terminal transport 
detection of a topological phase transition and Major ana 
zero- modes, along with braiding of Majoranas to harness 
their non-Abelian statistics?*™*' Such benefits provide 
strong motivation for pursuing Major ana physics in gate- 
defined wires in HgTe or related platforms. 

The remainder of this paper is structured as follows. 
Section [TT] explores the physical properties of the HgTe 
wires. We then turn in Sec. Ill to the application of 



Majorana zero-modes in this setting, treating both the 
clean and disordered cases within a simplified framework. 
Section [IV] summarizes our main results and discusses 
future directions in greater detail. Finally, three appen- 
dices contain additional calculations that further support 
the claims in this paper. 



II. CHARACTERIZATION OF GATE-DEFINED 
WIRES 

In this section we perform a detailed characterization 
of gate-defined HgTe wires. Section II A develops an 
analytical description of the system, starting from the 
Hamiltonian for a 2D HgTe quantum well and then sys- 
tematically including the effects of a confinement poten- 
tial, Rashba coupling, and applied magnetic fields. More 



accurate numerical simulations are explored in Sec. II B 
Our objective below is to demonstrate that such wires ex- 
hibit large and exceptionally tunable Rashba spin-orbit 
coupling and ^-factors as claimed in the introduction, 
rendering them promising for applications that will be 



briefly discussed in Sees. Ill and IV 



A. Analytic treatment of the confinement problem 

We begin by reviewing the well-studied physics of uni- 
form quasi-2D HgTe quantum wells, following conven- 
tions used in Refs. [36] and |38| Excitations in the four 
bands closest to the Fermi level can be described with 
a spinor ^(r) T = [^+(r), ^#+(r), (r), $ H - (r)]. 
Here ^^(r) and ^h,s(y) annihilate states with oppo- 
site parity at position r = (x,y) in the quantum well 
plane. Under time-reversal T these operators transform 
as ^e/h,s s ^e/h-s- The standard model Hamilto- 
nian for the 2D quantum well reads^^ 



2D 



J \ h* J 



(1) 



with 



h = [-M2D + D(dl + d$)]I + d a 
d = [-iAd x , iAdy, M + B{d 2 x + dl) 



M[eV] 


-0.01 


A[eVA] 


3.65 


B[eVA 2 } 


-68.6 


D[eVA 2 } 


-51.2 


R/(eE z )[k 2 ] 


-1560 


9e 


22.7 


9h 


-1.21 



TABLE I: Parameters for a 2D HgTe quantum well with thick- 
ness d — 70A. 38 At such a thickness the system happens to 
realize a topological insulator, though this property is incon- 
sequential for the formation of gate-defined wires in the bulk 
of the quantum well. 



In the above equations / is a 2 x 2 identity matrix, a 
denotes a vector of Pauli matrices, fi2D represents the 
chemical potential, and A, £?, D, and M are materials 
parameters dependent on the quantum-well thickness d. 

For the remainder of this subsection we focus on long- 
wavelength physics where it suffices to set B = D = 
0, as doing so greatly facilitates analytic treatment of 
the problem. Our numerics in Sec. |IIB| restore these 
terms to experimentally relevant values and confirm that 
they do not change the physics qualitatively. With this 
simplification Eq. ([!]) describes massive Dirac fermions 
with band energies 



E±(k) = -fi 2D ± Vm 2 + (Aky 



(3) 



The gap 2M for HgTe wells is quite small — typically on 
the order of O.OleVP^l As emphasized in the introduc- 
tion the Dirac structure together with this small mass 
cause the properties of gate-induced confined states in 
the bulk of HgTe to differ dramatically from those in 
conventional semiconductors such as GaAs or In As. We 
comment further on such distinctions below. 

Suppose that one now couples the quantum well to 
a set of top and bottom gates as shown in Fig. [lja). 
[The essential physics exhibited by this system can also 
be captured in the simpler experimental setup of Fig. 
[TJb), which contains onl y a sin gle bottom gate; we dis- 



cuss this further in Sec. |II A 1 ] These gates allow one 



(2) 



to separately tune the global chemical potential, the per- 
pendicular electric field in each region,^ and the relative 
potential between inner and outer regions which will de- 
fine a quantum wire of width W. Throughout this paper 
we assume that the voltages on the left and right pairs of 
gates are tuned identically to fully deplete carriers from 
the outer regions of HgTe. For now we also assume that 
each pair of top and bottom gates is adjusted symmet- 
rically so that structural inversion (z — >• —z) symmetry 
is present. (This restriction will be relaxed below when 
we discuss Rashba coupling.) If V(y) denotes the con- 
finement potential defining the wire, then under these 
conditions the Hamiltonian becomes H = B.2D + ^conf 



2 





| B-field | | 



FIG. 1: Schematic experimental setups, (a) A 2D HgTe quantum well of thickness d sandwiched by top and bottom gates and 
impinged by a magnetic field. The outer regions are tuned into an insulating state while the central gates define a quantum 
wire. Due to HgTe's unusual band structure, properties of the quantum wire (such as Rashba coupling and Zeeman splitting 
generated by a perpendicular magnetic field) are exceptionally tunable, (b) Simplified experimental setup featuring a single 
bottom gap. In this more e xperimentally feasible device, parameters for the wire exhibit essentially the same level of tunability 
as discussed in Sec. II A 1 (c) Replacing the central top gate by a proximity coupled s-wave superconductor allows the gate- 
defined quantum wire to host Majorana zero-modes. This application is especially attractive given the large spin-orbit energies 
(ru 30K) and giant effective ^-factors (which can reach ~ 600) that are achievable. 



where 



conf 



drV(y)* f *. 



(4) 



model the confinement with V(y) = 
\y\), though a more realistic smooth confine- 



Here we 
-VG(W/2- 

ment potential will be treated numerically later in Sec. 
II B| As illustrated in Fig. [2ja) states localized along y 
can exist in an energy window AE = min(|M|, \ V\). One 
can derive an effective ID Hamiltonian for these confined 
subbands — which we label by an index n = 1,2,.. . — by 
projecting the 2D quantum well Hamiltonian using 



#(r) -> 



[ z ikxX \ 

J k x \ 







$ n +(k x ,y) 




(5) 



The operators ^ n ±(k x ) above correspond to Kramer's 
pairs and annihilate states in the gate-defined wire with 
momentum k x in band n. One can obtain the two- 
component wavefunctions & n ± — (4>n,E±-,4>n,H±) T and 
associated band energies in the standard way by solving 
the Hamiltonian separately in the three regions of Fig. 
[2|a) and then matching boundary conditions (see Ap- 
pendix [A] for details) . Below we simply highlight some 
salient features of the problem. 

First, unlike for a conventional parabolic 2D dis- 
persion, the x- and ^-directions cannot be treated 
independently — hence the wavefunctions 3> n ,±(fec 5 y) de- 
pend on k x . We define overall phases such that $> nyS are 
purely real (which is always possible due to the form of 
H); moreover, these functions satisfy 



<f> n+ (k x ,y) = (-l) n a z <S> n+ (-k x ,-y) 
§ n+ (k x ,y) = $n-(-k x ,y). 



Note that except at k x = <$> n ,s(k x ,y) does not have 
well-defined parity under y — ^ —y. It follows from 
the properties above that inversion sends ip n ,s(kx) ~^ 
( — l) n iJj niS (—k x ) while under time reversal ^ n ,s(k x ) —> 
sl Pn,—s{~k x ) . 

Generally, increasing W reduces the energy difference 
between the confined bands, which scales like ~ ttA/W 
for \V\ ^> \M\. Increasing the depth V of the confining 
potential shifts these bands down in energy, allowing new 
confined states to emerge from the upper half of the Dirac 
cone. For V > 2M the lowest-energy confined bands 
begin to merge with the lower half of the Dirac cone; 
states in these bands remain confined at 'large' k x but 
are extended at 'small' k x due to hybridization with bulk 
states. This feature will be important in our numerics 



(6) 



discussed in Sec. |IIB| Within our analytical treatment, 
however, we avoid this complication for simplicity. 

Figure |2|b) illustrates the confined band energies and 
wavefunctions for a d = 70A quantum well hosting a gate- 
defined wire of width W = 800A and potential depth 
V = 0.025eV. To generate these curves we employed pa- 
rameters from Table |TJ which lists various quantities rel- 
evant for 70A-thick HgTe sheets^. These values will in 
fact frequently be adopted in our simulations below since 
for this thickness quantitative estimates are known for 
most parameters of interest to us here. As an important 
aside, we note that for the quoted ratio of M/B in the ta- 
ble HgTe resides in the topological insulator phased We 
stress, however, that whether the topological or trivial 
state appears has little bearing on the existence of con- 
fined gapless ID states that we seek to generate in the 
system's bulk. (Of course changing d to enter the trivial 
phase modifies the parameters in Table [I] and therefore 
has a quantitative effect on properties of the confined 
states. But the important point is that there is no sharp 
distinction in the two cases insofar as these levels is con- 
cerned.) 
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FIG. 2: (a) Schematic illustration of the gate- defined quan- 
tum wire. (b) Main plot: Band structure for a W = 
800A gate-defined quantum wire residing in a d = 70A HgTe 
quantum well. The materials parameters used appear in Table 
[I] while the depth of the confining potential is V = 0.025eV. 
Note that each band is doubly degenerate. Insets: Illustration 
of the wavefunction components <j) n ,E/H+ hi the n — 2 band. 



Upon expanding the confined band energies to second 
order in k x , one obtains a simple effective ID quantum 
wire Hamiltonian which serves as the starting point for 
our analysis in this paper: 



h 2 kl 

2m„ 



4>n, 



(7) 



where e n and m n respectively denote the subband en- 
ergy and effective mass for band n, \i is the effec- 
tive chemical potential, and the sum over pseudospin 
s = ± is left implicit. A rough estimate for the effec- 
tive masses of bands far from extended states can be 
obtained by setting k y = nn/W in Eq. p| and then 
expanding the Dirac spectrum to order /cj; this yields 
m n ~ (\M\/A 2 )^/l + [nirA/(MW)} 2 . Inserting param- 
eters from Table [i] and assuming W = 800A, one finds 
that |M|/A 2 w 0.006m e , where m e is the bare electron 
mass, while [kA/(MW)] 2 « 2. For these parameters the 
effective mass therefore increases appreciably with n. As 
an illustration, m n is ~ 0.01m e for n = 1 (which is com- 
parable to the effective mass for electron-doped InSb) but 



is enhanced to 0.04m e for n = 5 (which is close to the 
effective mass for electron-doped GaSb). Gate-defined 
wires in systems with conventional parabolic bands, by 
contrast, exhibit masses that to a first approximation are 
independent of the band index. 

Because of the symmetries imposed so far, at a 
given momentum each band is doubly degenerate. For 
applications — e.g., the pursuit of Majorana fermions — 
it is highly desirable to lift this degeneracy via pertur- 
bations such as Rashba spin-orbit coupling and applied 
magnetic fields We turn now to incorporating these 
ingredients into our effective ID Hamiltonian. 



1. Effective Rashba coupling 

Suppose now that the top and bottom gates in Fig. 
[TJa) are adjusted asymmetrically. Such asymmetric gat- 
ing generates a perpendicular electric field E z , resulting 
in a voltage drop U = E z d across the quantum well width 
d. The loss of structural inversion symmetry leads to 
Rashba spin-orbit coupling, which for simplicity we as- 
sume is induced uniformly throughout the 2D quantum 
well (this is by no means essential). Our objective here is 
to explore the effective Rashba coupling felt by the con- 
fined electrons in our gate-defined wire. We will begin 
with the regime where the electric field E z is weak (in 
a sense to be quantified below) so that one can extract 
this effective coupling within first-order perturbation the- 
ory. This perturbative analysis provides rough order-of- 
magnitude estimates for the achievable spin-orbit ener- 
gies characterizing the wire. We emphasize, however, 
that for 'large' E z modifications of the confined wave- 
functions produce striking non-perturbative effects which 
underlie our main findings in this paper. An initial dis- 
cussion of non-perturbative effects is provided below; ad- 



ditional results appear in Sec. |IIB| and Appendix |C] 

Our perturbative analysis begins with the known 
Rashba Hamiltonian for the 2D HgTe quantum wellp^ 



/ 



cZ 2 r* f 







R{-d x + id v ) 



\ 





R(d x + id y ) 



V 











/ 



(8) 



Here R = FeU, with e > the magnitude of the electron 
charge and F a material (and geometry) dependent pa- 
rameter. Upon projecting Hr onto the confined bands 
using Eq. {5}, one obtains an effective Rashba Hamilto- 
nian for the gate-defined wire, 

Hr^Hid,r= / ^[-«rnn'(t)ri + i'-+^4 (9) 

Jkx nri 

which contains both intraband and interband couplings 



4 



of strength 



>(k x ) = R J dy(/> niE +(k x ,y)(k x -d y 



)(/>n' t E-(k x ,y). 

(10) 

Time-reversal symmetry requires r nn >(k x ) = —r n > n (—k x ) 
while properties of the wavefunctions in Eqs. ([6| dictate 



that r nn >(-k x ) = (-1) 



'(k x ). Thus intraband 



Rashba couplings must be odd in k x , as are interband 
couplings that mix bands with n and n' differing by an 
even integer; all other interband couplings are even in 
k x . For the moment we will assume that the electric 
field E z is sufficiently weak that all r nn > with n ^ n' 
are small on the scale of the confined subband separa- 
tion and can hence be ignored. Continuing to focus on 
long-wavelength, low-energy physics, we expand the re- 
maining intraband couplings as r nn (k x ) « ha n k x and 
neglect terms of order k% and higher. Within these ap- 
proximations, in the presence of Rashba coupling the ID 
wire Hamiltonian in Eq. ([7]) becomes 



/i 



h 2 k 2 x 

2m n 



ha n k x cF y I ^n-(H) 



The Rashba coefficients a n follow from Eq. (10) and 
take the form a n = c n R/h, where c n are generically 
order-one dimensionless constants. Together with the ef- 
fective masses, these parameters define a characteristic 
spin-orbit energy for band n via Eso,n — \m n a^ l . Re- 
calling that R = Fell one can express this energy scale 
in terms of the voltage drop U across the quantum well 
as 



Es0 ' n -2 mn {—h~) ■ 



(12) 



To obtain rough numerical estimates consider a d = 
70A quantum well for which F = w -22. 3A (us- 

ing E z = U/d and R/eE z as given in Table 0. For sub- 
bands with effective mass m n ~ 0.01m e , a voltage drop 
U ~ 0.05V then yields a characteristic spin-orbit energy 
of Eso,n ~ 10K. Such scales reflect a roughly order-of- 
magnitude enhancement compared with s pin-o rbit ener- 
gies in electron-doped InAs or InSb wires pMU 

Let us now quantify the range of U over which the per- 
turbative analysis above holds. The physics is more uni- 
versal for high subbands whose minimum is far in energy 
from the extended bulk states, so we focus on such cases 



for simplicity. By inspecting Eq. (10) one sees that in- 



terband couplings that mix adjacent subbands scale like 
r n ,n+i(k x = 0) ~ R/X yi where \ y is the characteris- 
tic wavelength along y. Roughly, X y corresponds to the 
Fermi wavelength for electrons in region II of Fig. [2ja) 
so that r nin +i(k x = 0) ~ RV/A = FeUV/A for large n. 
Since the subband spacing scales like ~ ttA/W, interband 
mixing is unimportant for voltage drops U satisfying 



\u\ < 



nA 2 



e\FV\W 



(perturbative regime). 



(13) 



For larger voltage drops interband mixing — not only 
with other confined bands but also typically with ex- 
tended states since M is rather small — becomes impor- 
tant. The result is a dramatic reshaping of the con- 
fined wavefunctions by the perpendicular electric field, 
which has surprising and potentially useful consequences. 
Specifically, upon increasing U away from the perturba- 
tive regime, the effective Rashba energy characterizing a 
given confined band does not monotonically increase as 
one might naively expect, but instead undergoes striking 
oscillations. 

The existence of these oscillations can be anticipated 
based on the following argument. Without Rashba cou- 
pling, the Dirac dispersion for the 2D quantum well along 
k y with k x = is sketched in Fig. [3^ a). For a given en- 
ergy there exists only a single pair of wavevectors ±k y , 
so that within the central region confined states are built 
from plane- waves e ±%kyV . Switching on Rashba coupling 
splits the 2D bands as in Fig. |3jb) and changes the sit- 
uation qualitatively. In particular, two distinct pairs of 
wavevectors ±k y i and ±k y 2 now yield the same energy — 
hence confined wavefunctions involve superpositions of 
two harmonics, e ±zkyiy and e ±zky2y . The difference in 
these wavevectors increases with U, i.e., k y \ — k y 2 oc U. 
Consequently, varying U changes the profile of the con- 
fined wavefunctions in an oscillatory fashion. This effect 
is visible in Figure [3jc), which displays the numerically 
computed probability amplitudes \(/>E/H±(k x = 0,?/)| 2 
versus U assuming parameters specified in the caption. 
Oscillations in the confined wavefunction tails — which 
are clearly seen in the figure — in turn produce oscillations 
in the effective spin-orbit energies (and other physical 
properties as we will see) characterizing the gate-defined 
wire. 

The above qualitative argument for the appearance of 
oscillations relies on the existence of two harmonics in- 
side the central region where the wire exists. In contrast, 
the effect of the Rashba coupling outside the wire is of 
minor relevance for the physical properties of the con- 
fined states. This circumstance allows one to simplify 
the experimental setup of our device with little physical 
consequence. In particular, replacing the bottom gates 
in Fig. [TJa) by a single gate as shown in Fig. [ljb) leaves 
independent control over the global chemical potential, 
the confinement depth V, and the voltage drop inside 
the wire. One merely sacrifices independent tunability of 
the Rashba coupling outside of the wire — which in any 
case is unimportant. 

Quantitatively capturing such effects clearly requires a 
more exact treatment of Rashba coupling. In Sec. |IIB| we 
expose the oscillations using exact numerical simulations 
of the gate-defined wire. There we show that the effective 
Rashba energy scale for confined subbands can be tuned 
from zero to a maximum of a few tens of Kelvin, and 
back down to zero with a moderate increase in U (see 
lower panels of Fig. [5| . This level of tunability is highly 
attractive for spintronics and other applications. 

To close this subsection we remark that the Rashba os- 
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(a) R=0 (b) R#) (c) |<Mk x =0)| 2 l</» H+ (k x =0)| 2 




FIG. 3: Band dispersion at k x = for the 2D quantum well Hamiltonian (a) without Rashba coupling and (b) with Rashba 
coupling. The important distinction between the two cases is that there are generically two distinct wavevectors yielding the 
same energy in (a) but four in (b). Upon gate-defining a wire, confined wavefunctions are then built from additional harmonics 
when Rashba coupling is present. The added harmonics result in oscillatory behavior of the confined wavefunctions that, 
in turn, produce oscillations in physical quantities for the gate-defined wire such as the effective ^-factors and characteristic 
Rashba spin-orbit energies, (c) Evolution of the probability amplitudes \cj) E / H ±(k x — 0,?/)| 2 for the n — 5 band as a function 
of the voltage drop U . Data correspond to a d — 70A quantum well hosting a wire of width W = 800A and confinement depth 
V = 0.06eV, and were obtained numerically as described in Sec. |IIB| The oscillations in the wavefunction tails shown in each 
plot underlie the oscillatory behavior of physical quantities noted above. 



cillations discussed above are not special to Dirac systems 
such as HgTe. Indeed, our qualitative explanation merely 
required the existence of multiple harmonics, which 
would arise even in a conventional parabolic dispersion 
with spin-orbit splitting. Related non-perturbative phe- 
nomena have been explored in conventional semiconduc- 
tors in Refs. [HUH 



quantum well, 3 



2. Effective Zeeman splitting 



Let us now introduce an applied magnetic field and 
investigate the effective Zeeman splitting imparted to 
the gate-defined wire. We focus throughout on mag- 
netic fields B = Bz directed perpendicular to the quan- 
tum well since this orientation yields the strongest ef- 
fect by far. The wire's effective Zeeman splitting de- 
rives from two physically distinct contributions — the or- 
dinary (spin-orbit-enhanced) Zeeman effect as well as a 
component due to the orbital part of the magnetic field. 
The latter is often justifiably neglected in treatments of 
wires. Here, however, we show that the orbital contribu- 
tion dominates as a consequence of the Dirac spectrum 
exhibited by HgTe. Initially we treat the case without 
Rashba coupling, but discuss the strong interplay be- 
tween Rashba effects and Zeeman splitting at the end 
of this subsection. 

Consider first the standard Zeeman term for the 2D 



j 



d 2 r& 



9e 

9h 










-9e 




*, (14) 



~9h J 



which contains ^-factors Qe/h f° r the E and H sectors. 
Equation ([5| once again allows us to project onto the 
confined bands of interest. Taking the weak- field limit 
where interband mixing is negligible and expanding to 
leading order in k X: we obtain 



Hz — » Hid z 



1? 



9u,z^bB 



(15) 



with 



9n,z = f dy<f>l + (k x = 0,y)( 9E ) $„+ (k x = 0, y) . 
J \ g H j 

(16) 

Following our usual approach we estimate this compo- 
nent of the effective ^-factor by considering a quantum 
well of thickness d = 70 A, for which g& = 22.7 and 
gn = —1.21 according to Table [TJ Given these values 
one generally expects g Uj z to be of order 10 — still en- 
hanced compared to the bare electron ^-fact or but much 
smaller than that for, say, an InSb wire j^MJ Fortunately, 
as noted earlier the usual Zeeman term constitutes a sub- 
dominant contribution to the wire's effective factor. 
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Since we are dealing with quasi- ID states possessing 
strong spin-orbit coupling, the orbital part of the mag- 
netic field masquerades as an effective Zeeman splitting 
for the confined bands. To incorporate orbital effects we 
select Landau gauge for the vector potential and replace 
d x —> d x — ieBy/h in the 2D quantum well Hamiltonian in 
Eq. 0. (We continue to set B = D = for simplicity, so 
only the A term is affected by this substitution.) The or- 
bital contribution to the wire's effective ^-factor can be 
similarly obtained by projecting the orbital terms onto 
the confined bands using Eq. ([5|, neglecting interband 
couplings as appropriate for weak fields, and expanding 
to leading order in k x . This procedure yields the follow- 
ing term in our effective ID Hamiltonian, 



H 



where 



9n,orb 



lD,orb 



AAm e 



E 



(17) 



h 2 



J dyy$l + (k x = 0,y)a x <f> n+ (k x = 0,y). 

(18) 

Appendix [B] demonstrates that the above integral can 
be performed exactly, yielding 



dyy$Z i+ (0,y)<T x $ ni +(0,y) 



2M 



(19) 



which, remarkably, is independent of the band index, wire 
width W, and confinement depth V . It is illuminating to 
express the final result in terms of the two-dimensional 
effective mass for carriers in a uniform HgTe quantum 
well, = h 2 \M\/A 2 . Upon dropping the irrelevant 

band index label we obtain 
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Notice that as M goes to zero, orbital effects produce a 
divergent effective ^-factor; this is reminiscent of the di- 
vergent diamagnetic response for gapless Dirac systems 
such as graphene (see, e.g., Ref. |48|) . According to Ta- 
ble |lj a d = 70A thick quantum well is characterized 
by a very light 2D effective mass ra^ « 0.0057m e — in 
turn leading to a giant effective ^-factor g or b « 350 that 
greatly exceeds the Zeeman contribution as claimed. 

Several comments are in order, (i) References [37] and 
[381 previously emphasized the importance of orbital ef- 
fects on the ^-factor for quasi-lD states in HgTe quan- 
tum wells, but in the context of quantum spin Hall edge 
states. The physics in the two cases is similar but not 
identical. In particular, bulk inversion asymmetry terms 
(which are absent in our treatment) are essential for the 
effect in the quantum spin Hall case^H (ii) The Dirac 
dispersion, strong spin-orbit coupling, and small gap pro- 
vide the key ingredients underlying the giant ^-factor 
captured above. It is the Dirac structure that allows for 
a finite correction linear in B, which can be seen from Eq. 
(18) together with the symmetry properties of the Dirac 



wavefunctions in Eq. ([6| . Spin-orbit coupling intrinsic to 
the 2D quantum well Hamiltonian ensures that this lin- 
ear correction lifts Kramer's degeneracy, and the small 
gap guarantees that this happens very efficiently. (Hi) 
It is instructive to contrast our results for HgTe with 
the behavior for gate-defined wires in systems exhibit- 
ing conventional parabolic dispersion, e.g., an electron- 
doped GaAs quantum well. There the analogue of Eq. 
(18) would vanish by symmetry if the confinement po- 
tential is symmetric under y — > —y, so that the leading 
perturbative orbital effect appears at second order in B. 
A linear term at k x ^ could still arise if the wire forms 
from asymmetric confinement, but such a term would not 
manifest as an effective ^-factor for the confined bands 
since spin degeneracy would remain unbroken (at least in 
the absence of spin-orbit coupling). 

So far in our discussion of effective Zeeman splitting 
we have entirely neglected Rashba spin-orbit interactions 
induced by a voltage drop U across the quantum well 
width. The results above still hold in the perturba- 
tive limit where Rashba coupling is weak, modulo small 
corrections coming from orbital effects induced by the 
Rashba terms. In the non-pert urbative regime, however, 
the interplay between Rashba and orbital magnetic field 
effects produces still more striking physics. As described 



in Sec. II A 1 large voltage drops generate order-one mod- 
ifications of the confined wavefunctions that, crucially, 
are oscillatory in U. The effective ^-factor for the wire 
arises from projecting the orbital magnetic field terms 
using these modified wavefunctions, and hence inherits 
their oscillations. This effect is by no means small as we 
demonstrate numerically in the following section and an- 
alytically in a simplified model in Appendix [C| In fact, 
as we will see below the effective ^-factor for the confined 
subbands is exquisitely tunable and can be adjusted by 
factors of several hundred (and modulated in sign) by 
varying U over moderate voltage ranges. 



B. Numerical results 

Next we complement our analytic treatment above 
with more accurate numerical simulations of the gate- 
defined HgTe wire. The purpose of these numerics is 
twofold. First, in the preceding subsections several sim- 
plifying assumptions were made to facilitate analytical 
progress — e.g., focusing on the long- wavelength limit, 
considering a step-like confinement potential, etc. Here 
we simulate the full 2D quantum well Hamiltonian H = 
B.2D + H con f + Hr + Hz [the respective terms are defined 
in Eqs. ([I]), Q, (J8|, and ( [l4| )] with these assumptions re- 
laxed. We continue to focus on a d = 70A thick quantum 
well hosting a gate-defined wire of width W 800A, but 
now for the first time include the B and D terms in B.2D 
using the parameter values quoted in Table [TJ A more 
experimentally realistic confinement potential V(y), in 
which the confinement walls are broadened over a dis- 
tance of ~ 200A, will also now be taken; for a sketch see 
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FIG. 4: (a) Smeared confinement potential used to define a 
wire in our numerical simulations, (b)-(d) Calculated band 
dispersions in various cases assuming a d = 70A thick quan- 
tum well, wire width W = 800 A, and confinement potential 
depth V = 0.052eV. In (b) U = B = so that both Rashba 
coupling and a magnetic field are absent. A voltage drop 
U = 0.05V is present in (c) leading to Rashba splitting of the 
confined bands. A magnetic field separates the Rashba-split 
bands, as (d) illustrates for a perpendicular field of strength 
B — 0.1T. Dense bands at the top and bottom of these plots 
represent extended bulk states. 



Fig. (4^a). And finally, we properly account for orbital 
effects of the perpendicular magnetic field B = Bz by 
sending d x — » d x — ieBy/h in all terms in H (including 
the B,D terms and the Rashba Hamiltonian). The sec- 
ond and more important goal of our numerics is to quan- 
titatively capture the non-perturbative effects of Rashba 
coupling alluded to earlier — in particular the prodigious 
oscillations in the effective Rashba energy and ^-factor 
characterizing subbands in the gate-defined wire. 

It is convenient to treat the quantum well as continuous 
along x (to take advantage of translation symmetry) but 
discretized along y so that one can describe the system 
in terms of an effective tight-binding lattice model that 
is readily simulated. The discretization is achieved in the 
standard way. One first expresses terms involving d y in 
momentum space and then replaces k y — » sm(k y a)/a in 
the continuum Hamiltonian, where a is the discretized 
model's lattice spacing. Finally, a partial Fourier trans- 
form ^(k x ,k y ) = ^Ysj e ~ ikyVjl ®( k x,yj) results in an 
effective hopping problem on an TV-site chain with sites 
labeled by yj . In all simulations we take N = 500 and as- 
sume periodic boundary conditions along y to eliminate 
unwanted edge effects.^ The lattice spacing is adjusted 
such that the width of the system along y is aN = 10W; 
the exponential tails of the confined wavefunctions are 
then well-resolved numerically. 

Figures |4^b) -(d) present typical band structures calcu- 
lated within the above scheme, assuming a confinement 



depth V = 0.052eV. Case (b) corresponds to U = B = 
where neither Rashba coupling nor a magnetic field are 
present. Hence all bands there are doubly degenerate. 
Notice that the n = 1 and 2 subbands overlap with ex- 
tended bulk states at small k x . In (c) the voltage drop 
across the quantum well is set to U = 0.05V, resulting 
in a pronounced Rashba splitting of the n — 3 and 4 
subbands. A perpendicular magnetic field of strength 
B = 0.1T is present as well in (d) and produces a clear 
separation between these Rashba-split bands. 

In such simulations we quantify the Rashba spin-orbit 
energy for subband n by considering the B = limit and 
defining Eso,n = \^n^F y n- Here the Rashba parameter 
a n is deduced from the slope of the dispersion at k x = 
while kp,n is the Fermi wavevector for subband n when 
the chemical potential resides at the k x = crossing for 
that subband [e.g., \i ~ — 5meV for n = 3 and /i « 5 
meV for n = 4 in Fig. [4^c)]. This definition of Eso,n 
reproduces our previous expression ^m n a^ in the case 
of a simple quadratic dispersion, but is more appropri- 
ate when significant non-parabolicities arise as is often 
the case here. To extract the magnitude of the effec- 
tive ^-factor for subband n, denoted g n , we equate the 
magnetic-field-induced splitting of the confined bands at 
k x = with \g n \fiBt3- (One can infer the sign of g n ana- 
lytically as discussed below.) For concreteness we use a 
field strength B = 0.005T for this extraction throughout. 
Note that it is difficult to meaningfully compute Eso,n 
and g n for bands that intersect extended states at small 
k x {e.g., n = 1,2 subbands in Fig. [4|, so below we will 
not quote spin-orbit energies and effective ^-factors for 
such bands. 

Figure [5] illustrates the dependence of the effective g- 
factor (upper row) and spin-orbit energy (lower row) on 
the voltage drop U and confinement depth V for the first 
four confined subbands. White regions correspond to 
U, V values where a given confined subband either does 
not exist, or intersects bulk states at k x = so that 
the quantities of interest can not be determined as re- 
marked above. The dashed line in the figures roughly 
indicates the values of U, according to Eq. (13), where 



the crossover between perturbative and non-perturbative 
Rashba coupling regimes transpires. Consider first the 
perturbative limit. At U = the ^-factor ranges from 
~ 370 for the n = 1 band to ~ 470 for the n = 4 band. 
Such enormous values again arise because of orbital mag- 
netic field effects, with deviations from our previous an- 
alytical estimate arising primarily from the B and D 
terms. The attainable Rashba energies in this regime 
are of order 10K in agreement with estimates from Sec. 

Far more interesting is the non-perturbative limit 
where oscillations in both quantities are visible. The 
following key features are worth highlighting. The os- 
cillation frequency increases with the band index, which 
is why the variation with U over the range shown is rela- 
tively minor for the n = 1 band but becomes increasingly 
pronounced in the higher subbands. More importantly, 
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FIG. 5: Color-scale plot illustrating the remarkable oscillations in the effective g-f actor (upper row) and spin-orbit energy (lower 
row) for the four lowest confined subbands. The horizontal axis represents the depth V of the confinement potential used to 
define the wire, while the vertical axis corresponds to the potential drop U across the quantum well due to a perpendicular 
electric field. Data were extracted from numerical simulations of a d = 70A quantum well supporting a wire of width W = 800A, 
assuming a smooth confinement potential. For a given band, the ^-factor and spin-orbit energy are well-defined only if that 
band supports confined states at k x = 0. This condition cuts off the plots at small and large V. Dashed lines indicate the 
crossover between perturbative and non-perturbative Rashba coupling regimes according to Eq. (13). 



moderate changes in U effect giant modulations in both 
the effective ^-factor magnitudes and Rashba energies — 
the former varying from zero to more than 600, the latter 
between and ~ 30K. 71 Both quantities oscillate with 
roughly the same period but do so out of phase. That 
is, for a given subband the ^-factor amplitude reaches a 
maximum when the spin-orbit energy is minimized and 
vice versa. 

Appendix [C] treats a simplified model that allows one 
to analytically capture the main features of these oscilla- 
tions for high subbands. As the calculation is somewhat 
lengthy we will not comment on the details here but in- 
stead simply note two important conclusions. First, our 
analytical study reveals that for high subbands the oscil- 
lation period in U is approximately given by 



A/7 



4ttA 2 
eFVW' 



(21) 



Notice that the perturbative regime identified in Eq. ( 13 ) 
persists to roughly one quarter of a wavelength of the os- 
cillations. And second, Appendix [C] demonstrates that 
the lines in Fig. [5] at which \g n \ vanishes are associated 
with sign changes for g n . Thus the effective ^-factors for 
the confined subbands are highly tunable both in magni- 
tude and sign. 

It is important to address how the oscillations depend 
on the width W of the wire. On one hand Eq. (21 ) illus- 



trates that the oscillation period decreases with W. But 



on the other, a shortened period cuts off the quadratic 
rise of the spin-orbit energy with U in the perturbative 
regime [see Eq. (12)] at a reduced value of U. The net 
effect is that wider wires yield smaller attainable Rashba 
energies. Consequently, if one desires to maximize the 
effective spin-orbit coupling, narrow wires are generally 
advantageous. We have confirmed numerically, however, 
that the magnitude of the ^-factor oscillations remains 
roughly constant upon increasing W — at least up to val- 
ues W « 1500 A. This is perhaps not too surprising since 
in the perturbative regime Eq. ( |20| ) shows that the dom- 
inant orbital contribution to the ^-factor is largely in- 
sensitive to both U and W (unlike Rashba coupling). 
Measuring the giant ^-factor oscillations experimentally, 
with magnitudes peaking at ~ 600, should thus be even 
easier in wider wires. 



In summary, gate-defined HgTe wires possess the fasci- 
nating property that their subband-dependent ^-factors 
and spin-orbit energies can both be tuned continuously 
over enormous ranges simply by changing the gate volt- 
age. However, since these parameters vary out of phase, 
one cannot maximize both simultaneously. Values of U 
leading to a 'compromise' where both quantities remain 
large are still possible, however, and we shall exploit such 
cases in the next section when discussing one particularly 
appealing potential application — the pursuit of Majorana 
fermions. 



9 



III. MAJORANA ZERO-MODES IN 
GATE-DEFINED HgTe WIRES 

When a wire with an odd number of channels acquires 
a bulk Cooper-pairing gap, the system can form a topo- 
logical superconducting state supporting protected Ma- 
jor ana zero- modes at its endpoints. The physics is inti- 
mately related to that of the Kitaev chain introduced in 
Ref.|49| A particularly powerful means of fashioning such 
a setup experimentally was proposed by Lutchyn et alW* 
and Oreg et a/P^ (see Refs. [50-54 for multichannel gen- 
eralizations). These authors showed that a topological 
phase can be engineered in spin-orbit-coupled wires that 
are subjected to a magnetic field and proximity coupled 
to an ordinary s- wave superconductor. The magnetic 
field opens up chemical potential windows where an odd 
number of channels are occupied as desired. Spin-orbit 
coupling meanwhile causes the spin to depend nontriv- 
ially on momentum in each partially occupied band — 
allowing an s-wave order parameter to open a full pairing 
gap even in such odd-channel regimes. 

Ideally, wires featuring both large ^-factors and spin- 
orbit energies are desirable for this proposal. The for- 
mer permits one to operate at relatively weak mag- 
netic fields — hence disturbing the parent superconductor 
weakly — while the latter (among other benefits) leads to 
enhanced robustness against disorder as discussed below. 
It is interesting to explore the formation of Majoranas in 
gate-defined HgTe wires since they offer the possibility of 
satisfying both criteria simultaneously. A superconduct- 
ing proximity effect can be induced in the HgTe wire 
using a setup similar to Fig. [TJc) in which the central 
top gate is replaced by an s-wave superconductor. Such 
a configuration does not allow independent tuning of the 
the electron density and Rashba coupling for the wire, 
though it is conceivable that one can enhance the tun- 
ability by, say, employing additional top gates adjacent to 
the superconductor. In any case we assume in our anal- 
ysis below that Rashba coupling strengths and densities 
similar to those captured in the previous section can be 
realized here as well. We will first treat the clean case 
by studying numerically the full 2D quantum well Hamil- 
tonian with proximity-induced Cooper pairing, and then 
discuss disorder effects within a simplified effective ID 
Hamiltonian. 



A. Numerical phase diagram 

It is useful to first explore the rough phase diagram, 
and achievable gaps in the topological regimes, for the 
device in Fig. [TJc) assuming the clean limit. To this 
end we follow the methods outlined in Sec. Ill Bl to ex- 
actly diagonalize the 2D quantum well Hamiltonian H = 
H2D+H con f+HR+Hz+H/i in a perpendicular magnetic 
field B = Bz. The first four terms were previously sim- 



wire, Rashba coupling, and Zeeman splitting — including 
orbital magnetic field contributions where appropriate. 
The last term encodes (crudely) the Cooper pair- 
ing inherited from the neighboring superconductor. We 
model this term by pairing opposite pseudospins from the 
E and H sectors with pairing strengths A E and A#, 



J 



dr(A E ^E^E- + A h i/>h+iI>h- + H.c.) . (22) 



For simplicity, we have assumed spatially uniform A E / H 
above and neglected other symmetry- allowed pairing 
terms. Although we will simulate the full 2D Hamilto- 
nian, it is useful to note that one can ascertain the effect 
of proximity-induced pairing on the wire by projecting 
Eq. (22) onto the confined subbands. Such a projection 
produces both intra- and inter-band pairing terms whose 
magnitudes depend on the wavefunctions and the precise 
values of A e /h- ^t zero magnetic field time-reversal 
symmetry dictates that both A E and Ah are real but 
fixes neither their relative amplitudes nor signs. Absent 
detailed microscopic modeling which we will not attempt 
here, we set 



(23) 



ulated in Sec. II B and reflect the kinetic energy for the 
quantum well, smooth confinement potential defining the 



throughout to minimize the number of free parameters.^ 
All numerics discussed below were carried out for an in- 
finite W = 800A wire using parameters from Table |lj 

Figure [6](a) illustrates the chemical potential windows 
(blue regions) as a function of the confinement depth 
V where the wire possesses an odd number of channels 
as required for topological superconductivity. The data 
correspond to A = 0, B — 0.15T, and a voltage drop 
U = 0.05V. The chemical potential resides within the 
magnetic-field-induced gap in the n = 1 band at the left- 
most strip, the n = 2 band at the next strip over, etc. 
Oscillations in the band-dependent effective ^-factors ac- 
count for the varying width of these strips; recall Fig. 

El 

Upon turning on A, a topological phase supporting 
Major ana modes appears beyond a critical magnetic field 
in these odd-channel regimes. Figures |6|b) through (d) 
show B — /jl phase diagrams at constant V cuts [dashed 
lines in Fig. [6ja)] using A = 0.25meV. The shaded 
regions represent topological phases, the boundaries of 
which correspond to the fields that close the bulk gap. 
As in other wire setups^^, the minimum required field 
follows from E Zeeman « A, where E Zeeman = \\g n \n B B 
is the Zeeman energy for the topmost partially occupied 
band. (Roughly, this is the field required to overcome 
interband pairing.) Beneath each topological region we 
also list the effective ^-factor and spin-orbit energy char- 
acterizing the uppermost band. These quantities are en- 
couragingly large in all plots — for reference one may com- 
pare to electron-doped InSb wires for which g « 50 and 
Eso ~ IK. The small field scale at which the topological 
phase sets in [~ 25mT in (b)] is also noteworthy consid- 
ering the sizable 0.25meV pairing energy assumed. 
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FIG. 6: (a) Windows of chemical potential \i (blue regions) in which the gate-defined wire exhibits an odd number of channels, 
assuming a voltage drop U — 0.05V and perpendicular magnetic field B — 0.15T. The horizontal axis denotes the confinement 
depth V defining the wire. Inside of these windows, introducing proximity- induced superconductivity yields a topological phase 
supporting Majorana zero-modes at the wire's ends, (b)-(d) Phase diagrams in the B—ji plane along the constant- V cuts labeled 
by dashed lines in (a). Shaded areas denote topological phases, with the bulk gap (in units of the pairing energy A = 0.25meV) 
indicated by the color scale. For each topological region we also display the topmost partially occupied subband's effective 
^-factor and spin-orbit energy — both of which are quite large in all cases. 



The magnitude of the bulk gap in the topological phase 
is indicated by the color scale in Figs. [6jb)-(d). Near the 
phase boundaries with the trivial state, k x = excita- 
tions always determine the minimal gap. In the interior 
of the topological regimes, however, the gap is set by 
finite- k x excitations near one of the Fermi points. Some- 
what counterintuitively, the minimum excitation energy 
here need not be set by the topmost partially occupied 
band. Rather, in some cases 'background' confined sub- 
bands yield the minimum-energy gap. This indeed occurs 
in Figs. [6jc) and (d) and is responsible for the generally 
smaller gaps present there in comparison to Fig.[6jb). To 
illustrate the physics, we note that cut (c) corresponds 
to the band structure displayed in Fig. Bid) where two 
'background' subbands cross at E « — 4meV. Interband 
pairing becomes appreciable near that crossing and con- 
spires to reduce the gap somewhat. Such effects are likely 
non-generic but should be kept in mind. 

A still more favorable experimental situation appears 
in Fig. [7^a). The data here correspond to a larger 
voltage drop of U = 0.08V and a confinement depth 
V = 0.035eV, with the n = 1 and 2 subbands partially 
occupied. While the large effective ^-factor and spin- 
orbit energy for the uppermost confined band are com- 
parable to those in Fig.[6^c), the gap protecting the topo- 
logical phase is significantly larger and decays much more 
slowly with the magnetic field over the interval shown. 
The enhanced robustness follows simply because 'back- 
ground' confined subbands do not limit the gap here [in 
contrast to Fig. [6^c)]. 



B. Disorder effects 

Lastly we discuss crucial effects of disorder on the topo- 
logical phase. In the presence of time-reversal symmetry, 
Anderson's theorem dictates that random potential dis- 
order does not degrade the superconducting gap for s- 
wave-paired systems SIM] This theorem does not apply, 
however, to the topological phase since its formation re- 
quires a finite magnetic field. The severity of the gap's 
degradation by disorder depends on the degree to which 
time-reversal symmetry has been broken. A useful way to 
quantify this is via the ratio of the Zeeman and spin-orbit 
energies, Ezeem&n/ Eso, f° r the highest partially occupied 
confined band. Working in the limit -E/Zeeman /E SO < 1 
is highly advantageous since here spins at the Fermi level 
feel the effects of the field only marginally. In this sense 
time-reversal symmetry is weakly violated, imparting the 
system with greater immunity against disorder! 41 * 56 * 

At the minimum fields required to access the topologi- 
cal states shown in Figs. [6] and [7| the Zeeman energy falls 
in the range ^z eeman ~ 2 — 3K (roughly the size of A). 
This rough scale together with the spin-orbit energies 
listed in the figures suggest that it is indeed possible to 
stabilize Majoranas in the coveted spin-orbit-dominated 
regime Zeeman/ 'Eso <^ 1- (Since the topological phase 
requires Ez eem an ^ A, one can always trivially access this 
regime by making the proximity effect poor. The key 
point is that for HgTe wires this remains feasible even 
with generous values of A and gaps exceeding IK.) The 
parameters for Fig. [TFa) appear particularly promising 
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8000 lattice sites). One can then implement the random 
potential as 
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FIG. 7: (a) Phase diagram for a HgTe wire with voltage drop 
U = 80mV, confinement depth V = 35meV, and pairing 
strength A = 0.25meV. The color scale indicates the gap 
for the topological phase. The combination of a large gap, g- 
factor, and spin-orbit energy render this regime highly favor- 
able experimentally. In particular, the large spin-orbit energy 
implies reduced sensitivity to disorder, (b) Gap as a function 
of the ratio of the superconducting coherence length and the 
mean free path £/Z m f p , obtained by solving an effective disor- 
dered single-band ID model for the topological phase realized 
in (a). The data correspond to a magnetic field B = 0.2T and 
effective pairing amplitude A e ff = O.lmeV, with the chemical 
potential situated in the center of the Zeeman gap. The upper 
curve corresponds to HgTe parameters relevant for (a), while 
the lower curve corresponds to InSb wire parameters. The 
larger gap for the former stems from the greatly enhanced 
spin-orbit energy. 



due to the large attainable gap. 

Next we provide a rough illustration of the advantage 
afforded by the large spin-orbit energies found above by 
modeling a gate-defined HgTe wire by an effective single- 
band ID Hamiltonian with random potential disorder. 
We caution that our results here are only meant to ex- 
pose general trends. For one, disorder in the 2D quantum 
well will generate randomness in quantities aside from 
the local potential. The neglect of other subbands is 
also certainly crude since the most promising cases iden- 
tified above correspond to multi-channel situations. To 
mitigate the effects of 'background' subbands we will use 
parameters relevant for the topological phase in Fig. [7^a) 
since there additional bands are at least unimportant in 
the clean limit. 

With these caveats in mind, consider the ID Hamilto- 
nian H = H + disorder, where 



dxip 



mi 

2m 



cbA eff (^ + ^_ 



[i — ihaa y d x + ^-gfisBa 2 ) i/j 



H.c.) 



(24) 



describes the clean wire with proximity-induced pairing 
and ^disorder encodes the random potential. We choose 
parameters for Hq to reproduce quantities relevant for 
Fig. ga) at B = 0.2T (\ma 2 = 21.2K, ha = 2.2eVA, 
and g = 165, which is slightly reduced from the value 
quoted in Fig.[7^a) due to nonlinear effects). The Hamil- 
tonian is most easily simulated upon mapping the prob- 
lem onto a discretized lattice model (here we typically use 



^disorder = ^ Vxty\$x 



(25) 



with x now labeling discrete lattice sites. We choose 
the disorder potential V x to exhibit Gaussian white noise 
correlations with V x = and V X V X > = S XjX 'W 2 . The 
disorder strength W can be related to the mean-free 
path Z m f p = vyt via r _1 = ^ L W 2 a7V(ep), where a = 
6. 46 A which is the HgTe lattice constant and and 
N(ep) respectively denote the clean-system Fermi veloc- 
ity and density of states at the Fermi energy. 

We have performed simulations of H with various dis- 
order realizations in the case A e ^ = O.lmeV, B — 0.2T, 
and \i = 0, corresponding to the topological region of 
Fig. [TJa) with the chemical potential lying at the center 
of the Zeeman gap for the topmost bandP^ The upper 
curve of Fig.[7^b) displays the resulting disorder-averaged 
gap as a function of the ratio of the superconducting 
coherence length and the mean free path £/Z m f p , where 
£ = Trvph/Aefi. 41 Also shown for comparison are results 
obtained with InSb wire parameters — m = 0.015m e , 
Ha = 0.23eVA, g = 50 — using the same magnetic field, 
chemical potential, and pairing amplitude (note that the 
InSb lattice constant a = 6. 48 A is almost identical to 
that of HgTe). Due to the smaller ratio Ez eemari /Eso, f° r 
HgTe the gap takes on nearly the full value of A e ff in the 
clean case and remains substantially larger than that of 
InSb at comparable £/Z m fp- 75 This is consistent with the 
qualitative argument discussed above. Of course what 
matters ultimately is whether one can successfully induce 
a proximity effect and access the topological phase ex- 
perimentally in the first place. The results above (which 
again should only be taken as a rough guide) further 
suggest that pursuing Majorana fermions in gate-defined 
HgTe wires is a worthwhile endeavor. 



IV. DISCUSSION 

Gate-defined HgTe wires, for a number of reasons, offer 
great potential for the pursuit of applications requiring 
the manipulation of electronic spin degrees of freedom. 
The parent quantum wells can be fabricated with quite 
high mobility (at least up to 1.5 x 10 5 cm 2 V _1 s _1 )PZl. 
Confined subbands in the gate-defined wire can ex- 
hibit giant effective ^-factors (exceeding 600!) and large 
Rashba spin-orbit energies measuring tens of Kelvin. 
Even more striking is the exceptional tunability of these 
quantities evident in Fig. [5] — both can be altered from the 
large values quoted above through zero in an oscillatory 
fashion by moderate variations in gate voltages. 

It is useful to summarize the origin of these effects. 
The small gap and large intrinsic spin-orbit coupling for 
HgTe together cause the orbital part of the magnetic 
field to enormously enhance the effective ^-factors for the 
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wire. Such orbital contributions would typically be negli- 
gible in weakly spin-orbit coupled systems but dominate 
the Zeeman splitting here. As for the oscillatory behav- 
ior in the wire's effective Rashba coupling and ^-factors, 
these are rooted in non-perturbative modifications of the 
confined wavefunctions by a gate-induced perpendicular 
electric field. In principle wires patterned in more con- 
ventional electron-doped quantum wells (e.g., GaAs) can 
also undergo such oscillations, at least in the Rashba cou- 
pling, though likely with much smaller amplitude. Simi- 
lar physics is in fact implicitly present in the wide wires 
studied in Refs. I3T1J341 Armed with these insights, we 
suggest that gate-defined wires in hole-doped quantum 
wells may exhibit very similar physics to those in HgTe. 
This would be interesting to explore in greater detail in 
future work. 

As one enticing application, we explored the prospect 
of employing HgTe wires to stabilize Major ana zero- 
modes. The giant accessible ^-factors and large spin- 
orbit energies lead to a number of potential advantages, 
notably the ability to access a topological superconduct- 
ing state at quite small fields (as low as tens of milliTes- 
las) and with a sizable gap that exhibits reduced sen- 
sitivity to disorder. Here we wish to comment further 
on additional advantages offered by gate-defined wires, 
which apply not just to HgTe-based structures but to 
any suitable two-dimensional electron gas. In particular, 
the formation of wire networks appears to be relatively 
straightforward in this class of systems, requiring only 
additional patterning of gates on the quantum well. 

There are at least two interesting applications one can 
envision with such wire networks. The first is an im- 
proved detection scheme for the onset of the topologi- 
cal phase and accompanying Major ana zero- modes via 
transport. Consider, for instance, the multi-terminal 
setup shown in Fig. [8^a). Leads 1 and 2 in the figure 
allow one to inject current into the ends of the super- 
conducting part of the wire to search for the hallmark 
Majorana- mediated quantized zero- bias a nomalyp^"^ as 
has been done in recent experiments! 22 * 23 ^ Several authors 
have pointed out, however, that the closing of the bulk 
gap at the topological phase transition may be difficult 
to resolve in such a measurement, because the wave- 
functions for the gapless excitations may have very little 
weight near the ends of the superconductor (hence pro- 
ducing a weak transport signal) r 3 " 65 * Measuring trans- 
port from lead 3, which impinges on the bulk of the 
superconducting wire segment, should avoid this com- 
plication entirely and provide important complementary 
information about bulk physics. In particular, observ- 
ing a collapse and revival of the bulk gap coincident 
with the appearance of a stable zero-bias peak (even if 
not quantized) would provide extremely strong evidence 
for Majorana zero- modes. Such an experiment should 
also be able to distinguish 'accidental' zero-bias peaks 
driven by disorder^"^ or smooth confinement!^. A sec- 
ond, longer-term motivation of gate-defined networks in- 
volves braiding for the observation of non-Abelian statis- 




FIG. 8: (a) Multiterminal transport setup allowing for a re- 
fined detection of Majorana zero-modes. A current through 
leads 1 or 2 allows one to probe the hallmark zero-bias 
anomaly stemming from Majoranas. To readily resolve the 
gap closure (and revival) that accompanies the onset of the 
zero-bias peak, transport can also be measured from lead 3 
which probes the wire's bulk, (b) Braiding of Majorana modes 
may be performed using a network of superconducting wires 
that can be fashioned simply by deposition of additional gates 
on the quantum well. 



tics and quantum information applicat ions gee ^ e g ^ 
the setup in Fig. (8^b). There the keyboard of side gates 
should allow one to locally tune between topological and 
trivial regimes in a given part of the junction and hence 
transport Majorana zero-modes along the network. 

Various spintronics applications are also worth inves- 
tigating in HgTe-based wires. Spin qubits and spin tran- 
sistors are two natural candidates that warrant further 
exploration^"^. Finally, it would be quite interesting to 
perform a similar analysis of gate-defined quantum dots 
in HgTe quantum wells, which may inherit the remark- 
able features of the wires explored here. 
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Em = Ej . 



(A6) 



Appendix A: Solution to the confinement problem 

This Appendix provides details for the analytic solu- 
tion of the confined wavefunctions and energies in a gate- 
defined HgTe wire. As in Sec. II A we consider a Hamil- 
tonian H = H^d + H con f as defined in Eqs. ([I]) and 
set the parameters ji2D = B = D = in H2D for sim- 
plicity, and assume the step-like confinement potential of 
depth V illustrated in Fig. |2|a). The solution proceeds 
by projecting onto confined states using Eq. ([5| and then 
solving the Hamiltonian separately in regions I, II, and 
III labeled in Fig.^a) subject to the boundary condition 
that the wavefunctions are continuous. (The wavefunc- 
tions follow from a first-order differential equation when 
B = D = 0.) Since the Hamiltonian is block diagonal it 
suffices to focus on eigenstates & n ^+(k x , y) of the upper 
2x2 block; eigenstates of the lower 2x2 block are related 
by time-reversal symmetry. 

Consider first region II, where the solutions are de- 
scribed by plane waves oc e ±%kyV . Without loss of gener- 
ality we take V > so that confined states emerge from 
the upper half of the Dirac cone. The most general so- 
lution for the confined wavefunctions in region II then 
reads 



A(k x -\- iky) 



ce 



-M+^M 2 + A 2 (k 2 x + k 2 y ) 

-ik y y [ A{k x -ik y ) 

-M+ ^M 2 + A 2 (k 2 x + k 2 y ) 



with corresponding energies 



(Al) 



E II = -V+^M 2 + A 2 (k 2 x + k 2 ). (A2) 

In region I, solutions are evanescent waves oc e Ky with 
k > 0: 



($„,+)/ = a e K * 



A(k x + k) 
-M + <j^/M 2 + A 2 (k 2 x - k 2 ) 



Ei = ay / M 2 + A 2 (k 2 -K 2 ). 



(A3) 
(A4) 



Here a = ±1 represents the sign of the energy for a given 
confined state. Similarly, in region III we obtain a solu- 
tion with ft — ^ —ft, 



($n,+)//i = de 



A(k x ft) 



-M + a y / M 2 + A 2 {kl-K 2 ) 



(A5) 



Imposing continuity of the wavefunctions at the end- 
points of region II (i.e., at y = ±W/2) yields a set of four 
homogeneous equations for the four constants a, b, c, d 
appearing above. A non-trivial solution for these parame- 
ters exists provided the determinant of the corresponding 
matrix vanishes. This condition can be expressed as 



- 2 cos(k y W)Zj r Z_k y K + sm(kyW) x 

x [{kl - k 2 )Z\ + (k 2 x + k 2 y )Z 2 _ - 2k 2 x Z + Z. 

with 



(A7) 



-M - 



M 2 +A 2 {kl + k 2 y) 



-M + (Ty / M 2 + A 2 (k 2 x -K 2 ). (A8) 



As a further condition, the energies obtained in each re- 
gion must of course be equal; hence 



V + ^M 2 + A 2 (k 2 x + k 2 ) = a^M 2 +A 2 (k 2 x -K 2 ) . 

(A9) 



The two conditions in Eqs. (A7) and (A9) are sufficient 



to determine k y and ft, which depend both on k x and 
the band index n. These parameters can be obtained 
numerically, yielding confined-band energies of the form 



E n (k x 



V + J M 2 + A 2 [k 2 + ky t „(k x ) 2 



(A10) 



Furthermore, by requiring normalization of the wavefunc- 
tions the constants a, b, c, d can then also be determined 
uniquely up to an unimportant overall phase. 



Appendix B: Effective ^-factor due to orbital effects 



Section III A 21 discussed the orbital contribution to the 
gate-defined wire's effective ^-factor, which followed from 
the integral in Eq. (19). Here we show how one can eval- 
uate this integral analytically, yielding a result that is 
remarkably insensitive to details of the confined states. 
The first important step in the calculation is to observe 
that for k x = the upper and lower components of & n ,+ 
(respectively denoted by <j) n ,E+ and <p n ,H-\- ) ar e related 
to one another. Indeed, by inspecting Eq. (Al) one can 
show that in region II of Fig. [2ja) (i.e., for \y\ < W/2) 
the wavefunction components satisfy 



d y <j) ni E+(Q,y) 



dy<l>n t H+(0,y) = 



-Akl 



-M+JM 2 + A 2 k 2 



>n,H+(0,2/) ; 



(Bl) 



M+ JM 2 + A 2 k 2 



-</>n,E+(o,y) • 



(B2) 
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Similar relations hold in regions I and III (i.e., for \y\ > 
W/2): 



dy<l>n,E+(0,y) = 



An 2 



-M + (j\jM 2 — A 2 k? 
-M + cf\/M 2 — A 2 k? 



0n,if+(O,2/) , 

(B3) 

n ,£+(O,2/) . 

(B4) 



Consider next the normalization condition at k x = 0, 

/oo 
dy [4> n , E+ (p, yf + <Pu,h+(0, yf] . (B5) 
-co 

The right-hand side can in fact be recast into a form very 
similar to the integral in Eq. (19) that we are trying to 
evaluate. Specifically, upon inserting a trivial factor d y y 
(which equals unity) under the integral and then inte- 
grating parts, we obtain 

/CO 
d2/[0 n ,£+(O, y)yd y </) ni E+(0, y) 
-co 

+ (l>n,H+(0,y)ydy(l)n,H+(0,y)]. (B6) 

Finally, we break the right side up into separate integrals 
over regions I, II, and III and employ the relations from 
Eqs. (B1)-(B4) to eliminate the derivatives. After some 
algebra the three parts of integration can be reconciled, 
yielding 



4M 
~A~ 



d# n ,£+(0, 2/)# n ,tf+(0, y) , 



(B7) 



which immediately proves Eq. (19) 



Appendix C: Oscillations arising from 
non-perturbative effects of Rashba coupling 

In this final Appendix we provide a detailed account 
of the effects of Rashba spin-orbit coupling in the limit 
where a 'large' perpendicular electric field impinges on 
the quantum well. (The meaning of 'large' is clarified in 
Sec. II A 1| ) More precisely, our goal is to understand 
the dramatic oscillations in the gate-defined wire's effec- 
tive ^-factor and spin-orbit energy induced by varying 
the voltage drop U generated by the field (recall Fig. 
[5j. Ultimately these features reflect strong modifications 
of the confined-state wavefunctions by the perpendicu- 
lar electric field, and it is therefore essential that one 
treats Rashba coupling non-perturbatively here. Since 
this poses a nontrivial analytic task we will study a sim- 
plified model that facilitates progress yet still captures 
the essential physics. 

First, as in Appendix [A] we neglect terms quadratic 
in momenta in the 2D quantum well Hamiltonian 
and assume a step-like confinement potential V(y) = 



— VQ(W/2 — \y\) to define the HgTe wire. Second, we re- 
strict our considerations to deep confinement potentials 
and high subbands, where the oscillations are most pro- 
nounced as Fig. [5] illustrates. In other words, we assume 
V > M so that the Dirac cone inside the confined re- 
gion II in Fig. |2|a) is, roughly speaking, strongly shifted 
relative the cones in the surrounding regions I and III. 
The Dirac mass M in region II then negligibly impacts 
the confined states and can be safely ignored. Third, 
the Rashba coupling in regions I and III does not signif- 
icantly influence the oscillations we aim to describe, so 
for simplicity we will retain Rashba coupling only within 
region II (which rather naturally provides the dominant 
effect on the confined states). 

The full Hamiltonian we treat is then H = H^d + 
H con { + where the terms on the right are defined 
in Eqs. 0, Q, and ([8|; given the assumptions above 
we set B = D = 0, M MQ(\y\ - W/2), and replace 
R — »• R@(W/2 — \y\) in the Rashba term. As a final sim- 
plification we content ourselves with solving H above for 
confined wavefunctions with momentum k x = 0. This 
suffices for capturing directly the oscillations in the ef- 
fective ^-factor but not the effective Rashba energy scale 
characterizing the gate-defined wire, for which one would 
also need information about finite-/^ states. Neverthe- 
less, indirect arguments for Rashba oscillations can be 
made as described below. 

In the following, we proceed as in Appendix [A] and 
discuss the confined wavefunctions and energies at k x = 
by treating regions I, II, and III of Fig. [2ja) separately 
and then imposing proper boundary conditions. Note 
that the form of H guarantees that each component of the 
k x = wavefunctions has definite parity under y — ^ —y. 
More precisely, one can show that 

Vt>E+ (y) , <j>H+ (y) , $e- (y) , <j>H- [y)\ 
= X[(l)E+(-y),-^H+(-y),-^E-(-y)^H-(-y)] , (Ci) 

where solutions with A = +1 and —1 correspond to 
Kramer's pairs. (For notational simplicity, here and be- 
low we suppress the band index n on the wavefunctions; 
furthermore, all wavefunctions implicitly refer to k x = 0.) 
We therefore need only solve explicitly for the wavefunc- 
tions in regions II and III since the form in region I follows 
from Eq. flCl). 

We begin with region II, where the k x = wavefunc- 
tions are eigenstates of 



/ 



-V Ad v iRd y 



-Ad„ -V 



iRd„ 



-V Ad, 



-Ad v 



(C2) 



Diagonalizing this matrix using a plane-wave ansatz oc 
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e ±ik v y yields th e f our energies 



R + + R 2 



-V 



^//,2 = 1^1 o ^ 



77,3 



-£ 



77,2 • 



(C3) 



Let us assume that the confinement potential V is pos- 
itive so that the confined states emerge from the up- 
per half of the Dirac cone. In this case only the upper 
branches £77,1 and £77,2 are relevant so we discard the 
others. Most crucially, due to Rashba coupling there are 
now two distinct sets of momenta d=fc y i and dzfc^? which 
are related via 



R + V4A 2 + R 2 



%2 



-R + V4A 2 + # 2 : 



(C4) 



that yield the same energy. This fact reflects the usual 
Rashba splitting of bands [illustrated in Fig. [3jb) for 
M 7^ 0] and is intimately related to the appearance of 
oscillations. 



Consider next the corresponding confined wavefunc- 
tions, which generically consist of superpositions of plane 
waves with all four k y values above (since momentum is 
not conserved along y). Using Eq. (CI) with A = +1 
along with Eq. ( |C4| ), wavefunctions in region II take the 
form 



^% + {y)^ 

<P E -(y) 
\4>H-(y) ) 



= a+ 
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' (R + V4A 2 + R 2 ) cos(k yl y) 
2A sin(fc y iy) 
-i(R + V4A 2 + R 2 ) sm(k yl y) 
y 2iAcos(k y iy) 



\ 



1 (-R + VIA 2 + W) cos(k y2 y) 
2Asin(k y2 y) 
i(-R + V4A 2 + R 2 ) sm{k y2 y) 



V 



— 2iAcos(k y 2y) 



(C5) 



/ 



r 



where a + and b+ are constants and the added superscript 
on the wavefunction components indicates that A = +1. 
[For the Kramer's partner (/>e/ H ± with A = — 1 one sim- 
ply swaps cosines and sines, i.e., sin(- ••) —> —icos(- • • ), 
cos(- — )—> i sin(- • • ).] 

In region III the k x = wavefunctions are eigenstates 

of 



Hin 




with ^777 



M 



-Ad v 




Since there is no Rashba coupling here the so lutio ns can 
be immediately read off from Eqs. (A5) and (A6): 



1 <tis+(v) 



\ 



(y) 

\<t>H-(v) J 

with energies 



As before a 



—c\Ak 
Cxi-M + cry/M 2 



\ 



III 



A 2 n 2 ) 
—cI\Ak 
d x (-M + aVM 2 - A 2 k 2 ) 



(C7) 



Em = (j^M 2 - A 2 k 2 . (C8) 
±1 follows from the sign of the energy, 



the superscript A in the wavefunctions labels Kramer's 
partners, and c\ and d\ are constants. 

The wavefunctions must be continuous at y = W/2, 
which leads to a set of four homogeneous equations for a\, 
b\, c\, and d\. Again, for a non-trivial solution to exist 
the corresponding matrix must have zero determinant. 
After some algebra, this condition can be written as 



tan 



(fcyi + k y2 )W 



4A 2 + R 2 
2 V M 2 - A 2 k 2 



(C9) 



which holds for either A = +1 or —1. The energies in 
regions II and III must also match, yielding 



(R + a/4 A 2 + R 2 ) -V = a\/M 2 



(CIO) 



C9| ) and (CIO) are sufficien t to de- 
k y 2 and k. According to Eq. (C9) so- 



Equations (C4j), 
t ermine k y \ 

lutions for different bands correspond to the intersec- 
tion points of the different branches of the tangent with 
the ^-dependent right-hand side. Hence, by chang- 
ing the Rashba coupling R oc U, the sum k y \ + k y 2 
for band n can only vary in the interval k y \ + k y 2 G 
[2(n — l)7r/W, (2n — l)ir/W] (for positive energies) or 
kyi+k y 2 G [(2n — l)7r/W, 2rar/W] (for negative energies) . 
The actual variation with U is typically much smaller, 
and for the high bands of interest here one can approx- 
imate the sum as a constant, k y \ + k y 2 ~ 2fco- On the 
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^ 0.8 : 




0.05 0.1 
U[V] 



0.05 0.1 
U[V] 



0.15 




0.05 0.1 
U[V] 



0.05 0.1 
U[V] 



0.15 



FIG. 9: Oscillating behavior for components of one of the 
k x = wavefunctions evaluated at y = W/2, corresponding 
to the edge of the gate-defined wire. In the plots U is the volt- 
age drop induced by a perpendicular electric field, r]^j H+ = 

4> E/H+ (y = W/2), and V + /H _ = i<t> + E/H _{v = W/2). The 
data were obtained numerically by considering the n — 5 band 
in a HgTe system with d = 70A, W = 800A and V = 0.05eV. 



other hand, with Eqs. ( |C4| ), (C8) and (C10), one can 
show that the difference yields 



kyl ky2 



(V + E)R 
A 2 



VR 



(Cll) 



In the last step we used E <C V as appropriate for the 
deep confinement potentials we are considering. The two 
momenta are therefore well-approximated by 



V 

k v i Kko-^R, 



k y 2 ~ k 



V 

2A* 



R. 



(C12) 



Next, we consider the behavior of the wavefunctions 
at y = W/2 and define Ve/h+ = ^E/H+iy = W/2) and 
We/h- = i^E/H-iy ~ W/2) (the factor of i is inserted 
so that all 77's can be chosen real, which we henceforth 
assume is the case). By inspecting Eq. (C5) and its 



time-reversed partner one sees that the 77's for a given 
Kramer's pair are related via 



t 1e/h+ ~ ^e/h- 



^e/h- 



t 1e/h+ • 



(C13) 



Equation (C5) also implies that, schematically, for A 



+1 these components vary with ^1,2 according to 



(k y iW\ fk y2 W 



+ . /kyiW 

Vh+/ E - ~ c 3 sm I — — ) + c 4 sin 



2 

ky 2 W 



(C14) 



It then follows from Eq. (C12) that the arguments of the 



sines and cosines are given by 9 

oscillations in the wavefunction components as a function 



*^±%WR, leading to 



of R with the same period AR ■ 



8tvA z 

vw 



Note that - + 



t Ie/h+ 



and Ve/h- are ou ^ °^ P nase w ^h one another. This im- 
portant property is evident in Fig. [9j which illustrates 
these oscillations for the n = 5 band in a system with 
d 70A, W = 800A and V = 0.05eV; the data we re ob - 
tained by numerically solving Eqs. (C4), (C9) and (C10). 
Crucially, because the wavefunctions are continuous at 
y = W/2, the amplitudes of the exponential tails in re- 
gion III (and by symmetry region I) follow the same os- 
cillations. Thus Rashba coupling can, quite remarkably, 
be tuned such that for certain wavefunction components 
the exponential tail is completely suppressed. 

Let us now use these results to capture oscillations in 
the gate-defined wire's effective ^-factor. We focus on 
the orbital contribution to the ^-factor since as we saw in 
Sec. II A 2 this generally dominates over the conventional 
Zeeman terms. Of interest then is the energy splitting of 
k x = Kramer's pairs in a given band in response to the 
orbital part of the field. Upon sending d x — >• d x —ieyB/h 
in the 2D quantum well Hamiltonian (again with B = 
D = 0), first-order perturbation theory gives a splitting 



2eAB 



/>+*_ 



h 1 dyy [^+(y)^H+(y) 

(y)^_(y) - ^-(y)^E-(y)] 



H+ {y)4> + E+ {y) 

(C15) 



Note that this expression only involves wavefunctions 
with the same A — off-diagonal matrix elements vanish 
for symmetry reasons. Orbital effects from the Rashba 
term are also neglected for simplicity since they provide 
a small correction. 

Con tributions from 'small' y in the integral in Eq. 
(C15) are suppressed both by the factor y appearing 



in the integrand, as well as nodes in the wavefunctions 
present in the confined region (recall that we are treat- 
ing high subbands). Thus the integral is dominated by 
W/2. One can then restrict the range of in- 



values 



y\ 



> 



tegration to this regime and express Eq. (C15) in terms 



of 77's. We then obtain, using symmetry properties of the 
wavefunctions under y — >• — 



AE, 



orb 



8eAB 
h 



(Ve+Vh+ 



)/ 



)/ dy ye- 
w/2 



(ci6) 

The appearance of ^-factor oscillations is now manifest. 
When t]e/ h+ are both large while Ve/h- are both small 
(e.g., at U ~ 0.08V in Fig. |9| the former wavefunction 
components extend appreciably in space and hence pro- 
duce a large contribution to the effective ^-factor. The 
Ve/h- components, by contrast, are sharply suppressed 
away from the gate-defined wire and contribute negli- 
gibly. Similar results hold in the opposite limit where 
7 1e/h+ are sma ^ an d Ve/h- are l ar S e ( e -9-> at /7 = in 
Fig. [9]), though the sign of the splitting reverses. If, on 
the other hand, all components are of similar magnitude 
a cancellation results and the orbital contribution to the 



^-factor vanishes. Hence, the products of 77's in Eq. (C16) 



effectively double the frequency of oscillation compared 
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to that of the individual wavefunction components. Set- 
ting R = Fell, the period of (/-factor oscillations in terms 
of the voltage U is therefore given by 



Note that in contrast to Eq. (C15) this integral involves 



AU 



4irA 2 
eFVW 



(C17) 



Note that the 1/F-dependence of the periodicity is in 
accordance with Fig. |5j As another check, inserting V = 
0.05eV and the parameters used in Fig. [5j one obtains 
a period AU « 0.19V that is comparable to that seen 
in our more accurate numerics. Deviations arise mainly 
from the B and D terms in the quantum well Hamiltonian 
which have been neglected here for simplicity. 

Consider next the effective Rashba spin-orbit coupling 
induced in the gate-defined wire. Our analysis has so 
far focused on k x = 0, and extending the results to fi- 
nite k x is nontrivial. Nevertheless, terms in the Hamil- 
tonian involving d x , which generate Rashba-induced en- 
ergy corrections ex k x as k x — )• 0, can be treated relatively 
easily in the limit of small but finite k x . These correc- 
tions can be extracted from the k x = wavefunctions 
already determined. (Finite-fee corrections to the wave- 
functions contribute only at higher-order in k x .) Two 
such terms exist: one arising from the A term in the 2D 
quantum well Hamiltonian [Eq. Q], the other appear- 
ing directly in the Rashba Hamiltonian [Eq. In the 
non-perturbative Rashba coupling limit of interest, the 
splitting arising from the latter is limited by the Rashba 
coefficient R while that generated by the former is limited 
by A. Thus provided R< A, which roughly corresponds 
to the voltage regime considered here, the A term gener- 
ally dominates. 

Focusing on this contribution, the Rashba splitting for 
a given confined subband reads 



AE K = 2Ak x J dy [4>P + (y)4>H+(y) + <P + H + {v)<P~ E+ {y) 

- 4>t- {y)r H - (y) - (y)4>E- (y)} ■ (cis) 



coupling of Kramer's partners with opposite A. Since the 
integrand is highly oscillatory in the confined region due 
to the high subbands considered here, it is again a good 
approximation to evaluate the integral only in the region 



y\ > W/2. With the aid of Eq. (C13), one can then 



formulate Eq. (C18) in terms of Ve/h±^ yielding 



AE R n8Ak x (r)+ + r)+_ 



■Ve-Vh+'s 



W/2 



(C19) 

The above expression depends implicitly on R through 
the r] terms, and vanishes linearly at small R as one would 
expect from perturbation theory. To see this consider the 
products rj E ~ + r]H_ and rf^_rfjj + . According to Fig.pl at 
small R one of the z/'s in each product is roughly constant 
while the other increases linearly from zero with R; thus 
AE R ~ R as R -> 0. 



More interestingly, Eq. (C19) captures the oscillations 



in the effective Rashba coupling for the confined sub- 
bands at larger R. The physics is cl osely related to our 
discussion of the ^-factor below Eq. (C16); in particular 
the oscillations again derive from electric-field-induced 
modulations of the confined wavefunctions' exponential 
tails. Let us simply highlight two noteworthy points 
here. First, the oscillati on pe riod in U is the same as 
for the ^-factor — see Eq. (C17). Second, the Rashba and 



^-factor oscillations are o ut of phas e, which can be seen 
by comparing Eqs. (C16) and (C19) and inspecting Fig. 
pi Whenever either Ve/h+ or Ve/h- becomes zero, for 
example, the Rashba splitting vanishes while the ^-factor 
attains a large magnitude. This somewhat crude treat- 
ment of Rashba oscillations recovers the essential features 
seen in our more reliable numerics in Fig. [5] 
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Far away from the quantum wire (in the y direction) our 
system contains a narrow region with opposite field direc- 
tion such that the overall flux vanishes as required in a 
cylinder geometry. Such a field configuration has a negligi- 
ble effect, however, on the confined states. 
Note that having Eso,n — does not mean that effects 
of Rashba coupling on the energy spectrum for confined 
states is absent. Rather, this condition implies that the 
leading Rashba term is cubic in k x rather than linear. 
In contrast to the effective Rashba coupling and ^-factors 
for the wire, the effective pairing amplitudes do not os- 
cillate with U. The distinction arises because the y- 
integration of the projected pairing terms is dominated by 
\y\ < W/2; hence they are weakly sensitive to the oscilla- 
tory confined wavefunction tails. 

Note that if and Ah have similar magnitudes but op- 
posite signs, a very small proximity-induced gap may arise 
due to cancellations. Such a situation is non- generic, how- 
ever, and is not expected to be relevant in practice. 
A smaller value of A e ff is chosen here simply so that with 
InSb parameters the topological phase does not appear too 
close to the critical point at a field of B — 0.2T. 
It is tempting to use measured values for mean free paths 
for HgTe quantum wells and InSb wires to attempt a fur- 
ther comparison. We avoid this since the relevant quanti- 
ties in the presence of a superconductor are not known. 
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